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Abstract 

Background: Breast cancer stem cells (BCSCs) have been recognized as playing a major role in various aspects of breast 
cancer biology. To identify specific biomarkers of BCSCs, we have performed comparative proteomics of BCSC-enriched and 
mature cancer cell populations from the human breast cancer cell line (BCL), BrCA-IVlZ-Ol. 

/Wef/jotfe; ALDEFLUOR assay was used to sort BCSC-enriched (ALDH-n) and mature cancer (ALDH-) cell populations. Total 
proteins were extracted from both fractions and subjected to 2-Dimensional Difference In-Gel Electrophoresis (2-D DICE). 
Differentially-expressed spots were excised and proteins were gel-extracted, digested and identified using MALDI-TOF MS. 

Results: 2-D DICE identified poly(ADP-ribose) polymerase 1 (PARPl) as overexpressed in ALDH-H cells from BrCA-MZ-01. This 
observation was confirmed by western blot and extended to four additional human BCLs. ALDH-I- cells from BRCAl -mutated 
HCC1937, which had the highest level of PARPl overexpression, displayed resistance to olaparib, a specific PARPl inhibitor. 

Conclusion: An unbiased proteomic approach identified PARPl as upregulated in ALDH-f, BCSC-enriched cells from various 
human BCLs, which may contribute to clinical resistance to PARP inhibitors. 
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Introduction 

Breast cancer stem cells (BCSC), which are defined by their 
capacity of self-renewal and their ability to provide to non- 
tumorigenic, more differentiated cells, are increasingly recognized 
as major actors in breast cancer [1]. They are thought to be 
involved in critical aspects of breast cancer biology including 
carcinogenesis, metastasis, resistance to treatments, and thus 
clinical recurrence [2,3]. Accordingly, an effective targeting of 
BCSC is expected to significantly improve breast cancer outcome. 
As we previously demonstrated, flow cytometry-based detection of 
aldehyde dehydrogenase (ALDH) activity, can sort a subpopula- 
tion of cells from patient-derived xenografts and human breast 
cancer cell lines (BCLs) that have BCSC properties [4,5], including 



tumorsphere-forming capacity, serial passages in NOD/SCID 
mice and recapitulation of the cellular heterogeneity of the original 
cell lines. 

Proteomics-based profiling technologies may provide an unbi- 
ased approach to identify differentially expressed proteins between 
BCSC and non-BCSC populations, which may allow a more 
accurate definition of the BCSC subpopulation, a better under- 
standing of BCSC biology and eventually the identification of 
novel targets for BCSC-specific therapeutics. Using 2-dimensional 
Difference In-Gel Electrophoresis (2-D DICE), we have compared 
BCSC-enriched ALDEFLUOR-positive (ALDH-F) to ALDE- 
FLUOR-negative (ALDH-) BrCA-MZ-01 human BCL and 
found a significant increased expression of Poly(ADP-ribose) 
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polymerase 1 (PARPl), a promising target for a recently developed 
class of anticancer compounds. 



viability was evaluated by propidium iodide and used to define 
the sorting gates. 



Materials and Methods 

Cell lines and culture conditions 

BCLs were obtained from American Type Culture Collection 
(HCC1937, MDA-MB-436), from coUections developed in the 
laboratory of Dr S. Ethier, Karmanos Cancer Institute, Detroit, 
Michigan, USA (SUM 149 and SUM 159; obtained from Asterand; 
UKhttp: / / www.asterand.com/ Asterand/human_tissues/ 
hubrcelllines.htm), and from Dr VJ. Mobus, University of Ulm, 
Ulm, Germany (BrCa-MZ-01) [6]. BRCAl and BRCA2 gene 
status was checked for each cell line by fuU sequencing, as 
previously described [7]. The cell lines were grown using 
recommended culture conditions, as previously described [5]. 

Aldehyde dehydrogenase activity detection and cell 
sorting 

The ALDEFLUOR kit (StemCell Technologies) was used to 
isolate the population with high ALDH enzymatic activity using a 
FACStar PLUS (Becton Dickinson), as previously described [4,5] . 
Briefly, cells were incubated in a specific buffer containing ALDH 
substrate, while a negative control sample was obtained by 
identical incubation in presence of a specific 50 mmol/L of 
diethylaminobenzaldehyde, a specific ALDH inhibitor. Cell 



2D-DIGE experiments 

Proteins from ALDH+ and ALDH- BrCA-MZ-01 cell 
populations were extracted using an urea-based buffer (7.5 M 
urea, 0,01 ug/1 DTT, 2.5 M thiourea, 12. ,5 mM glycerol, 
62,5 mM Tris-HCl 2,5%) containing 1,5 mM of a protease 
inhibitor cocktail (Sigma-Aldrich, USA) and an equal amount 
(50 |J.g) of each sample was labeled with either Cy3- (ALDH—), 
Cy5- (ALDH+) or Cy2- (internal standard) CyDye according to 
the manufacturer's recommended protocols (GE Healthcare, 
Piscataway, NJ, USA), and pooled. Samples were then separated 
by 2-Drmensional gel electrophoresis with the following steps: 
isoelectric focusing, using 24-cm immobilized pH gradient strips 
(IPG 3-10 NL, GE Healthcare) and 2D-separation into 10% 
vmiform polyacrylamyde gels. Gels with CyDye-labeled proteins 
were digitalized using the Typhoon Trio Image Scanner (GE 
Healthcare), images were cropped with ImageQuant software (GE 
Healthcare) and further analyzed using the DeCyder v 6.5 
software package (GE healthcare). Spots of interest were excised, 
digested by trypsin and subjected to mass spectrometry analysis 
using MALDI-TOF MS (Ultraflex, Brucker Daltonics, Billerica, 
USA), using reflectron and positive modes with an ion acceleration 
of 25 keV. To process obtained mass spectra, we used the 
FlexAnalysis 2.0 software (Bruker Daltonics). To obtain protein 
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Figure 1. 2D-DIGE analysis of ALDH+ versus ALDH— BrCA-WlZ-OI cells. Upper panel. Representatives three-(left) and two-(rlght) color 
merged images of 2D-DIGE gel: red spots are from the Cy5-labeled ALDH+ cells, green spots are from Cy3-labeled ALDH— cells and blue spots are 
from standardized samples. Lower panel. Three-dimensional simulation of the relative abundance of PARPl protein in ALDH+ (a), standardized (b) 
and ALDH- samples (c). 
doi:1 0.1 371/journal.pone.01 04302.g001 



PLOS ONE I www.plosone.org 



2 



August 2014 | Volume 9 | Issue 8 | el 04302 



PARP1 and Breast Cancer Stem Cells 



A 



MDA-MB436 BrCA-MZ-01 SUM 149 SUM 159 HCC 1937 

ALD+ ALD- ALD+ ALD- ALD+ ALD- ALD+ ALD- ALD+ ALD- 




PARP-1 
Tubulin 



MDA-MB 436 BrCA-MZ-01 SUM 149 SUM 159 HCC 1937 

Figure 2. Up-regulated expression of PARP1 in ALDH+ relative to ALDH— cells from human breast cancer cell lines. A. Western 
blotting image. Presented blots are representative of at least 2 independent experiments. B. The quantitative comparison of PARP1 expression 
between ALDH+ and ALDH- cells, expressed as ratio. PARPl protein expression was first normalized to Tubulin expression. * PARP1 protein 
expression was compared between ALDH+ and ALDH— cells using Wilcoxon signed rank test. 
doi:1 0.1 371 /journal.pone.01 04302.g002 



identification, we used peptide mass fingerprint with an in-house 
Mascot server (Version 2.2.0, Matrix Science Inc., London, UK) 
probing the International Protein Index (IPI) human database 
fi-om the European Bioinformatics Institute, as described in [8]. 

Western blot experiments 

Protein lysates were loaded into SDS-PAGE, transferred on 
nitrocellulose membrane, blocked 1 h at room temperature in 
Tris-BufTered Saline/5% non-fat dry milk/0,1% Tween20, and 
incubated overnight with primary antibodies in blocking solution 
(PARPl and aTubulin mouse monoclonal antibodies, Sigma- 
Aldrich, USA). After extensive washings in TBS/0,1% Tween20, 
membranes were incubated 1 h at room temperature (RT) with a 
HRP-conjugated secondary antibody before being revealed with 
an enhanced chemUuminescence substrate (West Pico, Thermo 
Scientific, USA). 



Olaparib treatment 

Twenty-four hours after seeding, medium was changed and 
MDA-MB-436, BrCA-MZ-01, SUM 149, SUM 159 and 
HCC1937 cells were grown during 72 hours either in the presence 
of an inhibitor of PARP, Olaparib (AZD2281, Euromedex, 
France) at 10 |J,M (according to their respective IC50%) or in 
the corresponding concentrations of DMSO (control). ALDH+ 
and ALDH— cells from treated and control cell lines were sorted 
as described above and absolute cell numbers were counted using 
trypan blue. 

Results 

We used ALDEFLUOR-assay to sort ALDH+ and ALDH- 
cells from the BRCAl -mutated BrCA-MZ-01 human BCL and 
compared protein lysates obtained from these subpopulations by 
2D-DIGE (Figure 1). Among the proteins with differential 
expression (Table SI), PARPl was found as the most up-regulated 
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Figure 3. Effect of olaparib treatment on BCSCs from BCLs according to PARP1 expression. Human BCLs were seeded in culture flask and 
Incubated 72 hours in olaparib- or DMSO-containing medium. At each time and for each condition, percentage of ALDH+ and ALDH— cells were 
determined using ALDEFLUOR assay and absolute number of cells was counted using trypan blue. Results are presented after normalization by 
numbers of cells in DIVISO-treated conditions and Log-2 transformation. * Mean values of ALDH+ and ALDH— cells were compared using T-test. 
doi:1 0.1 371/journal.pone.01 04302.g003 



one in ALDH+ cells (ratio ALDH+/ALDH- = 1.56, 20 peptides 
matched, 23% of sequence coverage). 

PARPI protein expression was comparatively examined by 
western blot analysis in ALDH+ and ALDH— cell subpopulations 
from four additional human BCLs, including 3 other BRCAl- 
mutated BCLs (HCC1937, MDA-MB-436 and SUM149) and the 
BRCAl/2-wMtype SUM159 BCL (Figure 2). Overexpression of 
PARPI in BCSC-enriched ALDH+ population was confirmed in 
all samples, with a ratio ALDH+/ALDH— ranging from 1.59 in 
MDA-MB-436 to 4.99 in HCC1937 ceUs. In spite of a limited 
number of samples, this increase approached statistical signifi- 
cance (p = 0.06, Wilcoxon signed rank test). 

To determine the potential impact of PARPI overexpression on 
the sensitivity of BCSC-enriched cell subpopulation to PARP 
inhibitors, we exposed BCLs to 10 (Im olaparib during 72 hours 
and monitored the absolute number of ALDH+ and ALDH— 
cells. As shown in Figure 3, whereas olaparib activity was not 
significantiy different between ALDH+ and ALDH— cells from 
MDA-MB-436, BrCA-MZ-01, SUM 149 and SUM 159, the 
absolute number of ALDH+ cells in HCC1937, which displayed 
the highest level of PARPI overexpression, significandy increased, 
suggesting resistance to olaparib. In fart, there was a positive 
correlation between PARPI overexpression in ALDH+ cells and 
their relative resistance to olaparib (Figure SI). 

Discussion 

One of the main implication of BCSCs in breast cancer biology 
may be their potential role in therapeutic resistance [9]. Indeed, 
resistance to treatment, including chemotherapy and/or radio- 
therapy, was largely described in BCSCs and was thought to be 
related to various causes such as increased expression of adenosine 
triphosphate-binding cassette (ABC) transporters, resistance to 
apoptosis, lower proliferation and improved DNA repair ability 
[3] . By an unbiased proteomic-based approach, we have found an 
overexpression of PARPI in BCSC-enriched ALDH+ cell 
subpopulations from various BCLs. PARPI, the most abundant 
isoform of the PARP superfamily, is a chromatin-associated 



protein which participates in various biological functions such as 
cell proliferation, apoptosis, malignant transformation, transcrip- 
tional regulation and DNA repair [10]. Thus, PARPI plays a 
major role in base excision repair of DNA single-strand breaks. 
When DNA is damaged, PARPI recognizes the lesion and is able 
to bind it. Then, PARPI is activated and catalyzes poly(ADP- 
ribosylation) of various nuclear proteins. Ultimately, proteins 
modified by PARPI activity will aggregate and activate major 
components of DNA repair pathways. Notably, recent data, 
including our own [1 1], have demonstrated an adverse prognostic 
impact of PARPI overexpression in breast cancer [12]. 

In the absence of PARP 1 , DNA single-strand breaks accumu- 
late and lead to DNA double-strand breaks, which are not 
appropriately repaired if the BRCA patlnvay is deficient or 
dysfunctional, a phenomenon called synthetic lethality [13,14]. 
Consequendy, PARP inhibition was recently studied in various 
cancers, including breast cancer, with contrasted results. In 
BRCAl- or Si?C^2-mutated breast cancer, significant but 
transient and inconstant objective responses have been observed 
to olaparib in some studies [15,16], while other results were less 
conclusive [17]. Moreover, after promising results in a randomized 
phase II trial [17], iniparib (the actual anti-PARP activity of which 
remains discussed) in combination with chemotherapy failed to 
demonstrate any survival advantage in triple-negative metastatic 
breast cancer [18], a subset of tumors sharing some similarities 
with 5_RC/4 -mutated breast cancer. In the present study, we have 
observed olaparib resistance of ALDH+ cells from the BRCAl - 
mutated HCC1937 BCL, which was associated with the highest 
level of PARPI overexpression. A potential explanation to this 
observation could be that PARPI overexpression improves DNA- 
repair capacit)' of BCSC, and thus favors resistance' to DNA- 
damaging treatments, including olaparib. A similar phenomenon 
was suggested by a previous study in glioma cells showing that 
CD133H- cells, which are supposed to be enriched in stem cells, 
expressed higher levels of the DNA repair protein 06-methylgua- 
nine-DNA methyltransferase (MGMT) than CD133-, which 
correlated with resistance to radiation [19]. It can be argued that 
overexpression of a drug target usually results in drug sensitivity 
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rather than resistance. However, if this holds true with some 
targeted therapeutics such as trastuzumab or endocrine therapy, 
there are also other examples of drug resistance associated with 
overexpression of the target, such as BCR-ABL with imatinib in 
chronic myeloid leukemias or DHFR with methotrexate in 
sarcomas [20,21]. In these latter examples, and possibly with 
PARPl, overexpression of the target could increase the drug 
concentration required to effectively inhibit it. Another possible 
cause for increased resistance in ALDH+ cells could be 
overexpression of ABC transporters [22], which may function as 
drug-efllux proteins for various substrates, including olaparib [23]. 
However, by comparing gene expression profiles previously 
obtained from ALDH+ and ALDH— populations of the studied 
cell lines [5], we did not identify any significant difference between 
expression levels of ABC gene family (Figure S2). Identical results 
were found when restricting the comparison to the ALDH+ 
population of single HCC1937 cell line versus other cell lines, 
although such an analysis has an obvious limited statistical power 
(data not shown). 

An important limitation of this study relies upon the small 
number of cell hnes examined. Actually, only in HCC1937 cells, 
which display the highest differential expression of PARPl 
between ALDH+ and ALDH— subsets, was observed a significant 
increase in BCSC population under olaparib. Thus, PARPl 
overexpression may impede drug effect only above a certain 
threshold and this mechanism of resistance may only apply to 
some but not all breast cancers. This is consistent with clinical data 
showing exquisite antitumor activit)' for PARP inhibitors only in a 
subset of patients, even with BRCA-mutated tumors [15-17]. 

Clearly, the clinical relevance of our observations remains to be 
established and it is not known whether such a potential 
mechanism is actually operating in some BRCA-mutated breast 
cancer patients treated with olaparib. However, our results suggest 
that some intrinsic properties of BCSC such as increased DNA- 
repair capacit)', may contribute to resistance to PARP inhibitors, 
providing a basis for the relative failure of their clinical 
development in breast cancer. 
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Supporting Information 

Figure SI E£fect of olaparib treatment on BCSCs from 
BCLs according to PARPl expression. ALDH-H/ALDH- 
ratios of PARPl protein expression were plotted against the 
number of ALDH-I- cells after olaparib treatment (normalized by 
DMSO-treated cells) in various BCLs. * Pearson correlation 
coefficient. 
(TIF) 

Figure S2 mRNA expression of ABC family members in 
ALDH+ versus ALDH— cells. Gene expression data were 
obtained using Affymetrix U133 Plus 2.0 human oligonucleotide 
microarrays as described in [5]. Left panel. Heat map showing 
gene expression of ABC family members in ALDII-I- (blai k box) 
versus ALDH— cells (white box). * indicates ABC genes previously 
described as involved in drug-efilux [22]. Right panel, vulcano plot 
showing log2 ratio of fold-change in gene expression of ALDH-I- 
versus ALDH— (x axis) versus - log 10 ratio of p-values of Mann- 
Whitney U (y axis). 
(TIF) 

Table SI Differentially expressed proteins in ALDH+ 
compared to ALDH— BrGA-MZ-01 cells, identified by 
2D-DIGE and MS. 
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